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ABSTRACT 

Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from 
ultraviolet (UV) to sub-millimeter (submm) wavelengths, we estimate the total emission from dust 
and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model- 
independent as possible. We collect and exploit these data in order to measure from the spectral 
energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by 
dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more 
sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for 
more accurate estimates of dust temperatures and masses. 

The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three 
decades (from 10 -2 2 to 10 5 ). The inclusion of SPIRE data shows that estimates based on data not 
reaching these far-IR wavelengths are biased low by 17% on average. We find that the dust/stellar 
flux ratio varies with morphology and total infrared (IR) luminosity, with dwarf galaxies having faint 
luminosities, spirals having relatively high dust/stellar ratios and IR luminosities, and some early- 
types having low dust /stellar ratios. We also find that dust /stellar flux ratios are related to gas-phase 
metallicity (log(/d U st//*) = — 0.66±0.08 and — 0.22±0.12 for metal-poor and intermediate-metallicity 
galaxies, respectively), while the dust/stellar mass ratios are less so (differing by « 0.2 dex); the more 
metal-rich galaxies span a much wider range of the flux ratios. In addition, the substantial scatter 
between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. 
Comparing the dust/stellar flux ratios and dust temperatures, we also show that early- types tend to 
have slightly warmer temperatures (by up to 5 K) than spiral galaxies, which may be due to more 
intense interstellar radiation fields, or possibly to different dust grain compositions. Finally, we show 
that early- types and early-type spirals have a strong correlation between the dust/stellar flux ratio 
and specific star formation rate, which suggests that the relatively bright far-IR emission of some of 
these galaxies is due to ongoing (if limited) star formation as well as to the radiation field from older 
stars, which is heating the dust grains. 

Subject headings: galaxies: general - infrared: galaxies - galaxies: ISM - dust, extinction - galaxies: 
evolution 
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1. INTRODUCTION 

Nearly half of the bolometric luminosity of the Uni- 
verse is channeled through the mid- and far-infrared (IR) 
emission of galaxies (e.g., Hauser & Dwek 2001). This 
spectral region also probes important physical properties 
of galaxies, such as their metal content, dust content, and 
cold gas content (e.g., Draine et al. 2007). To understand 
galaxies we must first understand the physical processes 
that regulate their evolution, including the formation of 
stars and the interstellar radiation field (ISRF), and the 
return of radiant energy from these stars into the inter- 
stellar medium (ISM). Herschel stands poised to make 
major breakthroughs in these areas, by mapping galax- 
ies in the far-IR with unprecedented spatial resolution. 

Here we perform an empirical analysis of the galaxies in 
the KINGFISH project (Key Insights on Nearby Galax- 
ies: a Far-Infrared Survey with Herschel; P.I.: R. C. Ken- 
nicutt) , an imaging and spectroscopic survey of 61 nearby 
(d < 30 Mpc) galaxies, chosen to cover the full range of 
integrated properties and local ISM environments found 
in the nearby Universe. KINGFISH is closely follow- 
ing the observing strategy of SINGS (Spitzer Infrared 
Nearby Galaxies Survey; Kennicutt et al. 2003), by map- 
ping galaxies in their entirety with PACS (Photodctector 
Array Camera and Spectrometer; Poglitsch et al. 2010) 
at 70, 100, and 160/im and SPIRE (Spectral and Pho- 
tometric Imaging REceiver; Griffin et al. 2010) at 250, 
350, and 500^im. 

The physical processes contributing to the energetic 
output of a galaxy can be probed by constructing its 
spectral energy distribution (SED). Most of the stellar 
light is emitted in the ultraviolet (UV) to near-IR do- 
main, with the short-lived, massive stars dominating the 
UV and the more numerous older stars the near-IR. Dust, 
produced by the aggregations of metals injected into the 
ISM by massive stars through stellar winds and super- 
novae, absorbs the stellar light and re-emits it in the IR 
and submillimeter (submm) domains. 

Our goal is to determine empirically how much 
starlight escapes galaxies and how much is intercepted 
by dust, as a function of other galaxy properties. In par- 
ticular, we measure how the ratio of dust/stellar flux is 
correlated with properties such as morphology, metallic- 
ity, total IR luminosity, dust/stellar mass, dust temper- 
ature, and star formation rate. The KINGFISH sample 
is ideal for such an analysis, because these nearby galax- 
ies have been extensively studied with a variety of tele- 
scopes and models. Many of the galaxies' properties are 
already well-determined, and others are now better con- 
strained with data from Herschel. A secondary goal is 
to determine which subset (if any) of galaxies observed 
by KINGFISH could be plausible local counterparts of 
galaxies dominating the extragalactic background light 
(EBL). 

While our approach is an empirical one, there are alter- 
native theoretical approaches as well. A variety of differ- 
ent models have been applied to galaxy SEDs (e.g., Fioc 
& Rocca-Volmerange 1997; Silva et al. 1998; Devriendt 
et al. 1999; see review by Walcher et al. 2011), and there 
has been much recent work attempting to model both 
the stellar and dust SEDs of galaxies, over a wide range 
of wavelengths (e.g., Johnson et al. 2007; da Cunha et 
al. 2008; Noll et al. 2009; Buat et al. 2011; Popescu et 



al. 2011). Some of the results of these studies can be 
compared to ours (see Section [4J. 

This paper is organized as follows. In the next section, 
we describe the KINGFISH sample, and how we measure 
the galaxies' flux densities, from UV to submm wave- 
lengths. We explain how we integrate the SEDs to obtain 
dust/stellar flux ratios in Section [3] Then in Section [U 
we present our results, showing how the dust /stellar flux 
ratios are correlated with various galaxy properties. In 
Section we discuss how some galaxies in the KING- 
FISH sample could be considered local counterparts of 
galaxies contributing to the EBL. We end with a sum- 
mary of our results and conclusions in Section [5] 

2. DATA 

Our sample consists of 62 nearby galaxies, of which 61 
are in KINGFISH (Kennicutt, Calzetti et al. in prep.), 
and the other is M 33, observed by the Herschel M 33 
Extended Survey (HERM33ES; Kramer et al. 2010). The 
sample is chosen to cover a large range of galaxy proper- 
ties, such as morphological type, luminosity, metallicity, 
star formation rate, surface brightness, extinction, gas 
mass, dust content, radiation field strength, and ISM 
environment (see Kennicutt et al. 2003); however, the 
sample is not complete with respect to these properties. 

2.1. Flux Densities 

18 galaxies in our sample are in the Spitzer Local Vol- 
ume Legacy (LVL) survej{3, and for these we use the 
global flux densities measured by Dale et al. (2009). 57 
of the galaxies are in SINGS, and for those that are not 
in LVL, we use the flux densities measured by Dale et al. 
(2007), or when applicable, the updated values presented 
by Munoz-Mateos et al. (2009b). 

We obtained UV data from the Galaxy Evolution Ex- 
plorer (GALEX; Martin et al. 2005; 1528 and 227lA 
wavelengths); optical data from either the Sloan Digital 
Sky Survey (SDSS; York et al. 2000; ugriz bands) or 
Kitt Peak {BVRI bands); near-IR data from the Two 
Micron All Sky Survey (2MASS; Skrutskie et al. 2006; 
JHK bands); mid- and far-IR data from Spitzer's In- 
frared Array Camera (IRAC; Fazio et al. 2004; 3.6, 4.5, 
5.8, 8 fim) and Multiband Imaging Photometer (MIPS; 
Rieke et al. 2004; 24, 70, 160 /xm), and 100 /Ltm from 
the Infrared Astronomical Satellite (IRAS; Soifer et al. 
1987), when available; and submm data for one third of 
the galaxies from the Submillimeter Common- User Bolo- 
metric Array (SCUBA; Holland et al. 1999; 450 and 850 

For the two galaxies that were neither in SINGS nor 
LVL, IC 342 and NGC 2146, we obtained flux densities 
from various sources. We obtained UV magnitudes from 
GALEX (using the values quoted in Gil de Paz et al. 
2007 for NGC 2146), and converted them to fluxes using 
the calibration in Morrissey et al. (2007). We obtained 
optical fluxes from Buta & McCall (1999) and Marcum et 
al. (2001), respectively. We obtained 2MASS fluxes from 
Jarrett et al. (2003), using the calibration in Cohen et al. 
(2003). We measured Spitzer fluxes from their images, 
and we add the IRAS lOO^rn flux for NGC 2146. Lastly, 
M 33 (NGC 598) lacked optical fluxes in LVL, so we 
computed them from the data in Massey et al. (2006). 

18 http : //www . ast . cam. ac .uk/IoA/research/lvls 
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To these data we add Herschel far-infrared data from 
SPIRE (250, 350, and 500 fj,m). We have obtained 
SPIRE observations of 61 KINGFISH galaxies, includ- 
ing six galaxies (NGC 4254, NGC 4321, NGC 4536, 
NGC 4569, NGC 4579, and NGC 4725) that were ob- 
served as part of the Herschel Reference Survey (HRS; 
Boselli et al. 2010a); and we add the SPIRE observa- 
tions of M 33, from HERM33ES (Kramer et al. 2010). 
NGC 1404 and DDO 154 were observed by SPIRE but 
not detected, and their MIPS fluxes appear to be due 
to background sources (Dale et al. 2007), so we include 
these galaxies but regard their dust fluxes and masses as 
upper limits. 

The SPIRE (as well as PACS) flux densities will be 
shown and described in the KINGFISH photometry pa- 
per, Dale et al. (in prep.). For details about the other 
flux densities, see Dale et al. (2007, 2009). The galaxy 
distances are listed in Table [TJ and for most of the galax- 
ies, they are the same as those in Kennicutt et al. (in 
prep.). For a description of the distance indicators, and 
references for the distance measurements, see Kennicutt 
et al. (2003; in prep.). Some of the KINGFISH galaxies 
have been examined in detail in recent papers, such as 
NGC 1097 (Beirao et al. 2010; Sandstrom et al. 2010), 
NGC 3077 (Walter et al. 2011), NGC 6946 (Murphy et 
al. 2011b), and NGC 1291 (Hinz et al., in prep.). 

Elliptical apertures are used for the photometry and 
are chosen to approximately encompass all of the optical 
and infrared emission of a galaxy. Typically, this means 
that the 3.6 /im image was used to create the aperture, 
since 3.6 ^m is the bandpass within which Spitzer is most 
sensitive and the stellar disk is most spatially extended, 
although in a few cases the far-IR 160 /im disk is more ex- 
tended. The same aperture was used at all wavelengths. 
The global flux densities exclude foreground stars and 
background galaxies. We estimate that uncertainties in- 
volving the apertures may introduce up to 0.1 dex errors 
to the flux densities; however, these errors are quite small 
and do not significantly affect the dust/stellar flux ratios 
used for our analysis (see Section[3]). For details, we refer 
the reader to Dale et al. (2007 and 2009, Tables 1). 

Finally, for most of the KINGFISH galaxies that are 
also in SINGS, we measured Ha fluxes. For these, we 
used Ha images that were obtained as part of the SINGS 
ancillary program, either at the 2.1 m Kitt Peak Na- 
tional Observatory (KPNO) telescope or at the 1.5 m 
Cerro Tololo Inter-American Observatory (CTIO) tele- 
scope (Kennicutt et al. 2003). As described by Calzetti 
et al. (2007), exposure times were typically around 1800 s 
and standard reduction procedures were applied. There 
were 13 additional galaxies for which we did not have 
Ha fluxes but which were included in Kennicutt et al. 
(2008), and we used their fluxes for these galaxies, which 
are marked in Table [T] (with the superscript e ). For the 
galaxies whose Ha fluxes we could compare, our fluxes 
are slightly higher than those of Kennicutt et al., by 
« 0.17 dex, on average. For NGC 598, we obtained Ha 
data from Massey et al. (2006), from which we measured 
a flux density of 2.4 ± 0.2 x 10~ 10 ergcm~ 2 s _1 , which 
is slightly lower (by 0.15 dex) than the value given by 
Kennicutt et al. For NGC 2146, which was in neither 
sample, we use the Ha flux measured by Marcum et al. 
(2001). 



2.2. Morphological Classifications 

The morphological types and metallicities of galaxies 
play important roles in galaxy evolution (e.g., Calura 
et al. 2008; Fontanot et al. 2009), and they have been 
found to be correlated with the dust and stellar proper- 
ties of galaxies (e.g., Draine et al. 2007). Nonetheless, 
morphologies and metallicities are notoriously difficult 
to accurately determine without significant biases, and 
at fixed stellar mass, galaxies still have a fairly wide dis- 
tribution of morphologies (e.g., Bamford et al. 2009) and 
metallicities (e.g., Tremonti et al. 2004). In light of this, 
in Section [4] we will analyze the dust /stellar flux ratios 
of the KINGFISH galaxies as a function of other galaxy 
properties, but we will mark the galaxies by their mor- 
phological types (described in this section) and metal- 
licities (described in Section I2.3[) such that trends for 
late-types and early-types, and metal-poor and metal- 
rich galaxies, can be distinguished. Even if particular 
galaxies are misclassified or have ill-determined metal- 
licities, we expect the relative morphologies and metal- 
licities within the sample to be sufficiently accurate for 
statistical purposes. 

All but four of the KINGFISH galaxies were also in 
the SINGS sample (Kennicutt et al. 2003), in which 
galaxies were selected to span a wide range of RC3 (de 
Vaucouleurs et al. 1991) morphological types. Nonethe- 
less, many more data have been accumulated about the 
nearby galaxies since these classifications were made, and 
the classifications can now be done more accurately and 
homogeneously. 

Buta et al. (2010) have recently classified a subset of 
the objects in the Spitzer Survey of Stellar Structure in 
Galaxies (Sheth et al. 2010), using 3.6/im images with 
good spatial resolution. 29 of these objects are in the 
KINGFISH sample, and for these galaxies, we use the up- 
dated morphologies. In Table [1] we note the 17 galaxies 
which have been classified slightly differently than pre- 
viously. Some galaxies have been found to have slightly 
earlier types than previously; NGC 584 and NGC 855 — 
ellipticals which now have faintly detected disks — are ex- 
ceptions, as is NGC 1482, which now appears to have 
more of the structure of an Sa, rather than an SO. In 
any case, the previous RC3 classifications and Buta et 
al. (2010) classifications are generally consistent. In ad- 
dition, examining many of the same galaxies, Kendall et 
al. (2011) recently found that spiral structures are usu- 
ally similar in optical and IR images. We divide the 
KINGFISH sample into three types: there are 17 dwarf 
and irregular galaxies (Sd and later type), 32 spirals (Sa 
to Scd), and 10 early-types (E and SO). 

Galactic bars may also be an important property, be- 
ing related to gas concentration, star formation, and dust 
heating in the central regions of galaxies (Sheth et al. 
2005; Engelbracht et al. 2010). Approximately half of 
the KINGFISH galaxies have strong bars (SAB and SB), 
although because of the crudeness of the bar strength 
classifications, it is difficult to robustly determine the 
bar dependence of the galaxy properties. When signifi- 
cant, we quote the bar dependence of the galaxies' dust 
and stellar properties in the Section 01 the dependence is 
usually modest at most. 

2.3. Gas-Phase Metallicities 
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In order to quantify the metallicity dependence of cor- 
relations with the dust/stellar flux ratios, we use the oxy- 
gen abundances measured by Moustakas et al. (2010). 
In particular, we have chosen to use the abundances 
based on the theoretical strong-line calibration of Kob- 
ulnicky & Kewley (2004), even though it yields overes- 
timates for some galaxies. The empirical calibration of 
Pilyugin & Thuan (2005) is more accurate for typical L* 
galaxies, though it yields underestimates for star-forming 
metal-rich galaxies; more importantly, it was calibrated 
using only HII regions in spiral and irregular galaxies, 
and it may be dangerous to extrapolate beyond this 
regime. (For a discussion of the effects of using differ- 
ent metallicity calibrations, see e.g., Kewley & Ellison 
2008; Calura et al. 2009.) We use the 'characteristic' 
(globally-averaged) metallicities of the galaxies, even for 
those with metallicity gradients, such as NGC 5457. 

Engelbracht et al. (2005) found that there appears to 
be an oxygen abundance threshold at 12 + log (O/H) ~ 
8.2, such that galaxies below this threshold tend to have 
weak PAH emission; Draine et al. (2007) similarly found 
low PAH mass fractions in metal-poor galaxies. In or- 
der to determine how such metallicity transition may 
be related to the dust and stellar properties of galax- 
ies, we use only the relative abundances of Moustakas et 
al. (2010), and split the KINGFISH sample into thirds, 
consisting of 'metal-rich' (i.e., highest O/H metallicity), 
'intermediate' metallicity, and 'metal-poor' (i.e., lowest 
metallicity) galaxies. Our results are not significantly 
dependent on which metallicity calibration we use; the 
Pilyugin & Thuan (2005) calibration yields similar rela- 
tive abundances for most of the sample. For the galaxies 
lacking prominent emission lines, we estimate approxi- 
mate metallicities based on the -B-band L — Z relation 
(Moustakas et al. 2010; c/., Tremonti et al. 2004). Note 
that this relation has substantial scatter for faint galax- 
ies, and this is another motivation for focusing on the 
relative abundances. 

The absolute metallicity scales we use to split the sam- 
ple are 12 + log (O/H) = 8.88 and 9.08 (or 8.29 and 
8.42 using the Pilyugin & Thuan (2005) calibration), 
but these absolute abundances should be treated with 
caution. When we refer to 'metal-poor' and 'metal-rich' 
galaxies, these are meant to be relative to other KING- 
FISH galaxies. 

3. DUST/STELLAR FLUX RATIO 

3.1. Motivation 

Our goal is to estimate the emission from stars and 
dust in a way that is empirical and as model-independent 
as possible. This allows us to exploit the diversity of data 
that have been accumulated for these galaxies, and to 
compare to methods involving SED models (e.g., Draine 
et al. 2007; da Cunha et al. 2010). As described in de- 
tail below, we compute a dust/stellar flux ratio for each 
galaxy by integrating the SED from mid-IR to submm 
wavelengths and from UV to mid-IR wavelengths, and 
then taking the ratio. This is a quantity that can be 
physically interpreted as the amount of emission being re- 
processed by dust grains (mostly small and large grains) , 
relative to the unobscured emission from stars (especially 
young massive O and B stars, as well as intermediate- age 
AGB stars). 



Galaxy SEDs, and dust/stellar flux ratios in par- 
ticular, are related to other properties indicative of a 
galaxy's evolution, such as metallicity and morphology 
(e.g., Groves et al. 2008; Fontanot et al. 2009), which 
were discussed in Sections E21 and Section l2~3l Note that 
the dust /stellar flux is similar, but not equivalent, to the 
dust/stellar mass ratio, discussed in Section I4T1Z1 The 
dust/stellar flux ratio is also similar to the dust/FUV or 
IR/FUV ratio studied by many authors (e.g., Meurer et 
al. 1999; Kong et al. 2004; Johnson et al. 2007; Boquien 
et al. 2009; Wijesinghe et al. 2011), but these quanti- 
ties specifically measure the attenuation of UV photons, 
while the dust / stellar flux ratio also accounts for dust ab- 
sorbing optical photons from older stellar populations. A 
galaxy's dust/stellar flux is related to its specific star for- 
mation rate, as we will show in Section 14.41 and its star 
formation history. 

Lastly, note that galaxy geometry (i.e., inclination) 
and differential extinction will affect these ratios within 
galaxies to some degree (e.g., Jonsson et al. 2010). A 
few galaxies in the sample are highly inclined, such as 
NGC 4594 and NGC 4631, but as we will show in the 
next section, their dust/stellar flux ratios do not appear 
to be biased (with suppressed stellar emission). Simi- 
larly, Dale et al. (2007) examined the IR/UV luminosity 
ratio of SINGS galaxies and did not detect a significant 
trend with disk inclination. 

3.2. Procedure 

Our procedure is as follows. We begin by compiling the 
flux densities from the UV to far-infrared from Dale et al. 
(2007, 2009). For the ugriz-band optical SDSS data ob- 
tained for 11 of the 17 LVL galaxies in KINGFISH, and 
the J5Ui?/-band data obtained for the other galaxies, 
we first convert the apparent magnitudes to flux densi- 
tiesEH We then use the flux densities in units of Janskys, 
at wavelengths ranging from 0.15 to 850 /im. We have 
measured the SPIRE fluxes at 250, 350, and 500 /im our- 
selves, using the same apertures that had been used at 
shorter wavelengths. The SPIRE photometry and global 
flux densities will be presented in Dale et al. (in prep.). 

For galaxies missing data or detections at UV or 
submm wavelengths, we attempt to extrapolate the 
SEDs, in order to more consistently compare all of the 
galaxies in the sample. The extrapolations may be uncer- 
tain, but they yield more accurate dust/stellar flux ratios 
than neglecting these regions of the SEDs. Nonetheless, 
the effects of these extrapolations are usually extremely 
small, although as noted below, the UV extrapolations of 
a few galaxies have an effect of > 10%. We stress that the 
dust/stellar flux ratios are largely determined by the stel- 
lar and IR peaks of the SEDs, which are well determined 
for the 62 galaxies in our fiducial sample. 

For galaxies with SPIRE detections but lacking detec- 
tions with SCUBA at 850 /im, we perform a linear fit to 
the SPIRE log flux densities (at 250-500 /mi), and ex- 
trapolate to 850 /im, and give this flux an uncertainty 
of 1.5 times the uncertainty of the 500 /im flux. (The 

19 For the conversion from mag to Jy, we used the 
following conversion for the BV RI-haxi& magnitudes: 
http : //ssc . spitzer . caltech . edu/tools/magtojy/ref .html; 
and the following for the ugriz-hand magnitudes: 
http : //www. sdss . org/dr5/algorithms/f luxcal .html. 
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choice of 1.5 times the uncertainty is arbitrary; if we were 
to double this uncertainty, it would not significantly af- 
fect the final uncertainty estimated for the dust/stellar 
flux ratio.) We have verified that this is an accurate ex- 
trapolation for galaxies with both SPIRE and SCUBA 
detections, and in any case, it contributes a negligible ef- 
fect to the dust/stellar flux ratios: the ratios are affected 
by 0.002 dex on average and 0.03 dex (6%) at most. 

The extrapolation at the UV end is slightly more im- 
portant, because these wavelengths are closer to the stel- 
lar peak of the SEDs of some galaxies than 850 /im is to 
the far-IR peak. For the galaxies only lacking GALEX 
detections in the FUV, we extrapolate linearly from the 
NUV- and {/-bands (or NUV- and it-bands), and give 
this flux 1.5 times the uncertainty of the NUV flux den- 
sity. Three galaxies lack UV data altogether, and we 
examined each of these individually. The turnover of 
NGC 3077's SED is at longer wavelengths (« 800 nm), so 
we extrapolated linearly from the u- and (/-bands for this 
galaxy and accordingly gave these UV fluxes larger un- 
certainties. We assume that NGC 1377's SED resembles 
that of similar SO's in the sample (such as NGC 1266's 
in Fig. [!}, which turn over steeply. Using the mean 
and variance of FUV-NUV and NUV-B colors of these 
galaxies, we estimated the UV end of this galaxy's SED, 
and the uncertainty of the resulting dust/stellar flux ra- 
tio. Dwarf and irregular galaxies, including those in our 
sample, tend to have a shallower UV slope (Dale et al. 
2007). Using the mean and variance of FUV-NUV and 
NUV-B of these galaxies, we estimated the UV end of 
NGC 5408's SED, and the uncertainty of the resulting 
dust/stellar flux ratio. The extrapolations only out to 
the bluest FUV point have a very small effect on the 
dust/stellar flux ratios. The extrapolations through both 
UV bands are more significant; without them, one would 
overestimate the dust /stellar flux ratios of NGC 3077, 
NGC 1377, and NGC 5408 by 12%, 13%, and 50%, re- 
spectively. In any case, these factors are still relatively 
small on a log scale. 

In order to demarcate "dust" and "stellar" emission 
in the SEDs, the simplest approach would be to im- 
pose a strict wavelength cut at 5 /an, which is motivated 
by the fact that most stellar emission occurs at shorter 
wavelengths and dust mostly emits at longer wavelengths 
(e.g., Engelbracht et al. 2008). Nevertheless, a small frac- 
tion of stellar emission occurs in the mid-IR, and for some 
galaxies, and for early-types in particular, the emission 
from dust does not dominate until longer wavelengths. 
In addition, some galaxies have strong polycyclic aro- 
matic hydrocarbon (PAH) features that make it difficult 
to determine a strict demarcation of the stellar and dust 
SEDs. 

Motivated by this, we attempt to estimate the stellar 
contribution beyond 5 /im, as some authors have done 
(e.g., Draine et al. 2007; Munoz-Mateos et al. 2009a). 
For most galaxies, we fit a power-law to the stellar SED 
from the -RT-band to 4.5 /jm, and extrapolate to estimate 
this extra contribution (see Figure [1]) . Two exceptions 
are NGC 1377, for which we extrapolate from if-band to 
3.6 /im, and DDO 165, for which we extrapolate from K 
to 8 /im but exclude the 3.6 and 4.5 /im fluxes, as these 
appear to be slightly enhanced. This is not an ideal so- 
lution, because dust contamination may occur even at 
2 /im (Mentuch et al. 2010), but for most of the galaxies 
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Fig. 1. — Spectral energy distributions of NGC 1266 and 
NGC 1316. The dust/stellar flux ratios of the galaxies, with and 
without the SPIRE data included, are listed in the lower left of 
each panel. The dashed line indicates the estimated stellar con- 
tribution that is integrated at A > 4.5 /im. The open points at 
850 (im are extrapolated from the SPIRE flux densities. See text 
for details. 

in our sample, the Wien side of the dust SED domi- 
nates and its contribution increases with wavelength at 
A > 4.5 /im. The effect of accounting for the stellar con- 
tribution at these wavelengths has a very small effect on 
the dust/stellar flux ratios, and is significant only for the 
early-types (lowers their /dust//* by < 0.25 dex, and by 
0.5 dex for NGC 1404) and IC 342, but the contribution 
for this galaxy is very uncertain, because of its uncertain 
near- and mid-IR fluxes. We assume that the non-stellar 
contribution at A < 4.5 /im is minimal (see Meidt et al. 
2011). 

To estimate the stellar and dust emission, which we call 
/* and /dust, we integrate over the SED at A < 4.5/im 
and add the mid-IR stellar contribution for the former 
and integrate the SED at A > 4.5/im and subtract the 
same contribution for the latter. The area under the 
SED is computed directly, using Jy and Hz as the units. 
In particular, for each galaxy, we first linearly interpo- 
late log / (log v) over the range of wavelengths. Then we 
integrate over f(v) (i.e., not in log space). 

rX— 850/im , » rbQ^m , » 

/dust _ Jx=4.5tim aV J" ~ h.htirn ^ /star 

f rX— 4.5am -, , p50um , 1 >■ < 

J * Jx=0.15nm + h.^m 

where /star is the estimated stellar contribution at A > 
4.5/im, and is usually relatively small. Note that /dust 
is closely related to what some call total infrared lumi- 
nosity, which we discuss in Section 14.11 To estimate the 
uncertainties, we simply assume that the errors of the 
flux densities have a Gaussian distribution, sample from 
these distributions 10,000 times, and compute the vari- 
ance around the mean dust/stellar flux ratio. Because 
the stellar and far-IR SED peaks are usually well deter- 
mined for these galaxies, the uncertainties of /dust//* are 
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Fig. 2. — Distribution of dust/stellar emission ratios of the 
KINGFISH galaxies. The histograms show the distributions with- 
out (blue dashed line) and with (red solid line) the SPIRE fluxes 
included. 

usually small, even for the few galaxies with gaps in their 
SEDs, such as NGC 5408. 

3.3. Resulting Dust/Stellar Flux Ratios 

Two example SEDs are shown in Figure [1] NGC 1266 
and NGC 1316 are both SOs, and as we shall show later, 
they have similar masses. Nevertheless, they have very 
different dust/stellar flux ratios. In addition, they are 
among the galaxies that, before Herschel, did not have 
any far-IR data beyond 160 ^m. 

The distribution of dust/stellar emission for our fidu- 
cial sample of 62 galaxies is shown in Figure [5] For com- 
parison, we also show the /dust//* distribution without 
the SPIRE data (and without the 850 /im extrapola- 
tion). In other words, we compare to the distribution 
one would have obtained before Herschel, where the only 
far-IR data were from MIPS (and IRAS and SCUBA, for 
a few galaxies). 

For some galaxies, especially those with SEDs peaking 
at long wavelengths, Herschel contributes stronger con- 
straints on the ratio of dust/stellar emission, by tracing 
additional cold dust components not detected by Spitzer. 
Consequently, by adding a contribution in the far-IR, 
the SPIRE fluxes slightly increase the dust/stellar flux 
ratio of some galaxies in the sample. The mean ra- 
tio is log /d ust//* = -0.59 ± 0.07 without SPIRE and 
log /dust//* = -0.52 ± 0.07 with SPIRE, an increase of 
about 17%E3 (Note that the standard deviations about 
these means are « 0.54 dex.) NGC 1512 is the galaxy for 
which the inclusion of SPIRE data has the largest effect, 
increasing its dust/stellar flux ratio by 0.23 dex. The 



20 The mean ratios, without using logarithms, are /dust//* = 
0.51±0.08 without SPIRE and f duBt /f, = 0.55±0.09 with SPIRE, 
an increase of ss 8%. Since /dust//* spans three orders of magni- 
tude, however, we use the logs of the ratios throughout this paper. 



mean ratio varies with morphology as well, as shown in 
Table [2] and discussed in Section [JJ 

The three galaxies with the largest dust/stellar flux 
ratios are NGC 1482, NGC 2146, and NGC 2798, which 
are all starbursting early-type spirals. NGC 1482 also 
has evidence of superwind outflows (Hota & Saikia 2005), 
while NGC 2146 and NGC 2798 have peculiar morpholo- 
gies, possibly indicative of interactions or mergers (but 
see Greve et al. 2006 on NGC 2146). The two galaxies 
with the smallest ratios are NGC 1404 and NGC 584, 
which are both massive ellipticals (or E/S0) with very 
little FIR emission; NGC 1404 is also experiencing ram- 
pressure stripping as it falls through the Fornax cluster 
(Machacek et al. 2005). NGC 1404's / du st//* (~ 0.006) 
should be considered an upper limit, because its FIR and 
submm fluxes appear to be due to a background source 
(Dale et al., in prep.), and it is an indication of the most 
stellar-dominated a galaxy's SED can be. 

In Table [TJ we show the dust / stellar flux ratios of the 
KINGFISH galaxies. The near-IR morphologies, TIR lu- 
minosities, dust masses, stellar masses, and star forma- 
tion rates are also listed in the table, and are discussed 
below. 

4. RESULTS: CORRELATIONS WITH GALAXY 
PROPERTIES 

4.1. Total Infrared Luminosity 

4.1.1. Estimating Ltir, 

We begin by analyzing the total infrared (TIR) lumi- 
nosity of the galaxies in our sample. The TIR luminos- 
ity is a useful quantity because it can be directly inferred 
from the IR fluxes, and because it can be used as a proxy 
for the obscured star formation as well as the tempera- 
ture of dust grains (e.g., Dale & Helou 2002, Draine & 
Li 2007). We follow Draine & Li (2007), and use the 8, 
24, 70, 160/xm data from IRAC and MIPS to estimate 
the TIR luminosity: 

Ltir = 0.95(2/1^)7.9 + 1.15{vL v ) 2 4, + (vL u ) 7 i + (vL^ieo 

(2) 

We have also tested the Dale & Helou (2002) formula, 
which uses only the MIPS bands, and have obtained 
very similar results. (We have chosen not to include the 
SPIRE bands in Eqn. [2] because a calibration of Ltir 
with these bands has not yet been developed and tested.) 

Estimates of Ltir are designed to encompass all of the 
emission from PAH particles, very small grains, and large 
grains, the proportions of which depend on the starlight 
density distribution and the relative abundances of the 
grain populations (Draine & Li 2007; Compiegne et al. 
2011). Therefore, it is useful to compare Ltir and the 
dust/stellar flux ratio, in order to analyze the emission 
from dust grains vis-a-vis stellar emission. 

4.1.2. Results 

The correlation between Ltir and the dust /stellar flux 
ratio of the KINGFISH galaxies is shown in Figure As 
mentioned in Section [3j the total IR luminosity is very 
closely related to the quantity /dust- Plotting /dust X D 2 
versus /dust//* yields a result very similar to that shown 
in the figure. 

In Figure [3^, the galaxies are also labeled by their 
morphologies, such that we distinguish galaxies that are 
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KINGFISH Galaxy Properties 
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The columns are: galaxy name; morphological type, from Kennicutt et al. (2003); redshift-independent distance (see Kennicutt et al., 
in prep.); dust/stellar flux ratio (Eqn. [TJ; TIR luminosity, using Draine & Li (2007) calibration; dust temperature, estimated from 
modified blackbody fit to far-IR SED; dust mass, converted from dust temperatures using Li & Draine (2001) 500/^m mass absorption 
coefficient; stellar mass, using Zibetti et al. (2009) calibration; star formation rate (SFR), from FUV and TIR luminosities, and 
from Ha and 24 fim luminosities. The mass and SFR errors are underestimates: they include only formal errors from the fluxes, not 
systematic errors. a Morphology obtained from Buta et al. (2010), and was different than that listed in Kennicutt et al. (2003). b 
Mdust estimated from MIPS fluxes only, without longer wavelength fluxes. c Stellar mass-to-light ratios estimated from g — i and 
i — H colors; the others are estimated from B — V and V — H colors (see Zibetti et al. 2009). d Far-UV flux extrapolated from longer 
wavelengths. e Ha flux obtained from Kennicutt et al. (2008). f Ha flux obtained from Marcum et al. (2001). 
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classified as E and SO, Sa to Scd, and Sd and later- 
type. The spiral galaxies NGC 2146 and NGC 1097 have 
the largest TIR luminosities, while the dwarf galaxies 
DDO 154, DDO 165, and M81 Dwarf B have the faintest 
luminosities. The late-type spiral galaxies NGC 598 and 
NGC 7793 have almost exactly the same /dust//* and 
Ltir; as noted by Smith et al. (1984), these galaxies 
have similar photometric and kinematic properties, but 
different spiral arm structures. 

Many of the spirals tend to be found in the locus of 
/dust//* just below unity and L T ir > 10 43 ergs _1 in 
the figure. Most of the galaxies with small dust/stellar 
flux ratios and large TIR luminosities are early-types, 
although this may be partly due to a selection effect, 
as these galaxies were selected to be detectable in the 
IR (Kennicutt et al. 2003). In contrast, those with 
small /dust//* and small Ltir, are dwarfs and irregu- 
lars. Here and in terms of other galaxy properties, some 
of the earlier-type spirals follow the trends of the E 
and SO galaxies, while others have properties more sim- 
ilar to other spiral galaxies. The trend in Figure [3] is 
more complicated than a simple morphological distinc- 
tion, however, as some early-types have relatively large 
dust/stellar flux ratios (/dust//* > 1); as we will show 
later, these galaxies also tend to have larger specific star 
formation rates and dust temperatures. 

Note that Dale et al. (2009) have a somewhat sim- 
ilar plot (their Fig. 3), in which Ltir/Lb is plotted 
on the horizontal axis, rather than the dust/stellar flux 
ratio. Their sample is obtained from the LVL survey, 
which was designed to be nearly volume-limited. As a 
result, it is dominated by faint dwarf galaxies, and some 
of these fill part of the parameter space in the lower 
right of the figure, with low Ltir and slightly higher 
dust/stellar flux ratios than the dwarf galaxies in KING- 
FISH. In comparison, KINGFISH is incomplete and has 
relatively few dwarf galaxies, so the survey mostly covers 
the high- TIR part of the parameter space. Lastly, Soifer 
et al. (1989) also performed a similar analysis, plotting 
infrared/ visible flux versus FIR luminosity for a sample 
of IR-bright galaxies in IRAS. They obtained a single 
trend albeit with significant dispersion, but most of the 
dwarf galaxies in KINGFISH would likely not have met 
their selection criteria. 

We summarize the morphology dependence of TIR lu- 
minosity and other galaxy properties analyzed in this 
section in Table [2j which lists the means and standard 
deviations of these properties as a function of Hubble 
type. These could be useful as a local benchmark for 
comparisons with high-redshift studies, such as studies 
of submillimeter galaxies (e.g., Santini et al. 2010). 

As mentioned above, approximately half of the KING- 
FISH galaxies have strong bars, and it is possible that 
the presence of a bar affects emission by dust and stars 
in the central region of a galaxy. We find that the 
mean of the log dust/stellar flux ratio is —0.57 ±0.10 for 
weakly barred galaxies (no bar or SA) and —0.48 ± 0.09 
for strongly barred galaxies (SAB or SB). On the other 
hand, for the total infrared luminosity, the means are 
log Ltir = 42.9 ± 0.2 and 43.3 ± 0.2 for weakly and 
strongly barred galaxies, respectively. This is consis- 
tent with the observation that most IR-selected starburst 
galaxies are barred (e.g., Hunt & Malkan 1999). Never- 



theless, the bar fraction has also been observed to be 
higher among galaxies with higher optical luminosities, 
because these galaxies became dynamically cool and suf- 
ficiently massive to host bars earlier than fainter galaxies 
(Sheth et al. 2008). Thus one might expect barred galax- 
ies to have brighter stellar and infrared luminosities than 
weakly barred ones. Our data in fact bear this out: the 
mean stellar and dust luminosities, which we quantify as 
/»(47rD 2 ) and /dust (47rD 2 ) , are both higher in barred 
galaxies by ss 0.43 dex, the same factor as log Ltir. 
Therefore, it is simply the case that strongly barred 
galaxies in KINGFISH are more luminous than weakly 
barred ones, and this effect cancels in the dust/stellar 
flux ratios. 

In Figure [3)d, we show the metallicity dependence of 
the relation between the TIR luminosity and dust / stellar 
flux ratio. Firstly, it can be seen that the metallicity de- 
pendence and morphology dependence are similar: dwarf 
galaxies tend to be relatively metal-poor, spiral galax- 
ies tend to have intermediate metallicities, and early- 
type galaxies tend to be metal-rich. Nonetheless, it is 
not a one-to-one relationship: some early-type spirals 
with large Ltir are also metal-rich (such as NGC 1097), 
while some intcrmediatc-mctallicity galaxies are early- 
types with relatively large /dust//* (such as NGC 855). 

In general, from the morphological and metallicity de- 
pendence of Ltir and /dust//*, we can tentatively infer 
an evolutionary sequence from the figure, such that as a 
typical late-type galaxy grows and becomes more lumi- 
nous, it becomes more metal enriched and has a larger 
dust fraction; this is also accompanied by more emission 
by stars, a larger metallicity, and a growing stellar mass 
(cither by star formation or a merger), as the galaxy 
becomes an early-type. For the early-types, either stel- 
lar mass growth outweighs dust production, because the 
dense molecular clouds in which dust is typically accreted 
have dissolved, or a substantial amount of dust grains 
are ejected or destroyed, such as by supernovae shock- 
waves and thermal 'sputtering' (e.g., Draine & Salpeter 
1979; Dwek 1998; Pipino et al. 2011). This interpreta- 
tion of Figure [3] is merely speculative, however, and as- 
sumes that typical present-day late-types resemble past 
stages of present-day early-types; that is, it assumes that 
the growth of galaxy disks precedes that of bulges (e.g., 
Bournaud et al. 2009; Ceverino et al. 2010), although 
there is evidence that this assumption may be too sim- 
plistic (e.g., MacArthur et al. 2009; Bundy et al. 2010). 

In contrast, the dwarf and irregular galaxies have lower 
TIR luminosities (most have Ltir < 10 42 ' 5 erg s _1 ) and 
relatively little emission from dust (usually /dust//* < 
0.3). They appear to be a distinctly different galaxy 
population, exhibiting different properties than typical 
late-type spirals and inhabiting different environments 
(e.g., Leroy et al. 2008; Gavazzi et al. 2010). 

In any case, although general trends are apparent in 
Figure [31 there are plenty of exceptions and variation 
within the KINGFISH sample. Note that the sample was 
selected not to have strong AGN, defined as an AGN that 
dominates substantial portions of a galaxy's spectrum, so 
it is unlikely that AGN contribute much to this variation. 

4.2. Dust and Stellar Mass 

Next, in order to add to this picture of galaxy evolu- 
tion, we can use the infrared SEDs of the galaxies in our 



Dust and Stellar Emission of Galaxies 



9 



45 



44 



K 
J" 

O 



43 



42 



40 





1 i i i i | i i i i | i i i i | i i 
ngc ^ 

NGC 1097 site* _ 


(a) 

- 




♦ 


alA zn 

* >\ 8 S3 

♦ CO 01 — 

* 3 




* NGC 855 










A dwarfs - 




* ♦ ■ spirals 




J early-types - 

i , , , , i , , , , i , , , , i , r 



1.5 -1 -0.5 

lo g f dust/ f * 



0.5 



45 



44 



0) 

s 

an 
o 



43 



42 



41 



40 



(b) 



± 



O 



J_ 



I 1 1 1 1 I ' ' ' 1 I 
NGC 1097 3146$ 



* 

s NGC 855 



5* 1 



zo 

°M 
05 
CO® 



_L 



_L 



▲ metal-poor — 

intermediate 
• metal-rich - 



■1.5 -1 -0.5 

lo g f dust/ f * 







0.5 



Fig. 3. — TIR luminosity estimated using Draine & Li (2007) formula (using 8, 24, 70, 160/^m data) versus dust/stellar emission. Left plot: 
red circles, green squares, and blue triangles have morphology classifications of E and SO, Sa to Scd, and Sd and later type, respect ively . 
Right plot: red circles, green squares, and blue triangles have the highest, intermediate, and lowest oxygen abundances (see Section 12.31 1: 
open points indicate abundances estimated from the B-band luminosity- metallicity relation (Moustakas et al. 2010). 



TABLE 2 

Mean Galaxy Properties as a Function of Morphological Type 



log /dust//* 



log M duBt /M t 



log Ltir 
(log erg s" 1 ) 



Tdust 

(K) 



g SSFR Ha + 24iim 

(log M Q yr- 1 ) 



log SSFR FVV+Tm 
(log M Q yr- 1 ) 



All Galaxies 
Dwarfs 
Spirals 
Early-types 



-0.52 ±0.07 ±0.54 
-0.69 ±0.09 ±0.36 
-0.35 ±0.07 ±0.39 
-0.87 ±0.28 ±0.88 



-2.95 ±0.09 ±0.68 
-2.70 ±0.13 ±0.53 
-2.83 ±0.08 ±0.49 
-3.77 ±0.26 ±0.83 



43.12 ±0.14 ±1.08 
41.88 ±0.26 ±1.08 
43.74 ±0.09 ±0.51 
43.07 ±0.18 ±0.57 



27.0 ±0.6 ±4.3 
27.6 ±0.9 ±3.8 
25.8 ±0.6 ±3.5 
30.2 ±2.0 ±6.0 



-9.84 ±0.10 ±0.71 
-9.52 ± 0.12 ± 0.43 
-9.90 ±0.12 ±0.65 
-10.15 ±0.37 ±1.04 



-9.74 ±0.09 ±0.74 
-9.27 ±0.11 ±0.47 
-9.76 ± 0.09 ± 0.53 
-10.45 ±0.34 ±1.06 



KINGFISH galaxy properties: 
35 spira ls (S a to Scd), and 10 
14.11 14.31 and 14.41 respectively, 
galaxies with Ha fluxes from 



mean±bootstrap error±standard deviation. There are 17 dwarf and irregular galaxies (Sd and la ter) , 
early-types (E and SO), /dust//*, M dust /M,, L Tm , T d 

list? and SFR arc discussed in Sections [3l 14.21 
SSFR refers to specific star formation rate, SF R/ M,\ the quantities in these two columns include 
Kennicutt et al. (2008) and galaxies with extrapolated FUV fluxes (see Section [3j) . 



sample to estimate dust and stellar masses. These are 
more physical quantities than fluxes, but they require 
some model assumptions, and hence have additional sys- 
tematic uncertainties. 

4.2.1. Estimating Td us t and Mdust 

We estimate dust masses using dust temperatures, 
which are determined from simple fits to the FIR SED 
(e.g., Hildebrand 1983). We use the MIPS and SPIRE 
FIR and submm (70 — 500/im) flux densities and per- 
form fits using a single temperature blackbody modi- 
fied by an emissivity law proportional to A - ^. Our ap- 
proach is similar to that of Engelbracht et al. (2010) and 
Gordon et al. (2010), who also applied this method to 
galaxies with MIPS and SPIRE data. The use of a sin- 
gle temperature fit for whole galaxies essentially yields 
an average dust temperature; in practice, most galaxies 
have multiple dust components at different temperatures, 
such as in photo-dissociation regions and in the diffuse 
ISM, although one component might dominate. We as- 
sume j3 — 1.5, which is consistent with recent observa- 
tional constraints at far-IR and submm wavelengths (e.g., 



Dunne & Eales 2001; Paradis et al. 2009; Gordon et al. 
2010). In any case, the resulting dust temperatures are 
slightly dependent upon the assumed emissivity, and it is 
in a systematic way, such that (3 — 2 would yield slightly 
lower temperatures by a few degrees for the whole sam- 
ple (e.g., Bendo et al. 2003). We also do not attempt to 
account for the 'submm excess' inferred by some authors 
which may be due to a very cold dust component or to a 
wavelength dependent emissivity law (e.g., Galametz et 
al. 2009; Gordon et al. 2010). We infer the uncertainties 
of the temperatures with a Monte Carlo analysis that 
includes the flux errors. 

Most of the resulting temperatures range from 20 — 35 
K, with NGC 1512 having the coldest temperature (w 
21 K), and NGC 1377 and NGC 2146 having excep- 
tionally warm ones (« 43 K and 37 K, respectively). 
NGC 1404 yields an extremely cold temperature, but 
its FIR and submm flux densities are very small, and 
the 70 and 160 /an fluxes quoted in Dale et al. (2007) 
appear to be due to a background source, so we discard 
the galaxy's temperature as unrealistically low. For most 
other galaxies, a modified blackbody provides a good fit 
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to the FIR SED. The temperatures and uncertainties are 
listed in Table Q] 

Given Tdust of a galaxy, we estimate the dust mass with 
the following: 

duSt " ^x(^B x (T dust )) [6) 

where f\ is the flux density, D is the distance from the 
galaxy, B\ is the Planck function, which is 2ckT/X 4 in 
the Rayleigh-Jeans limit (and the additional 4ir factor is 
due to integrating over steradians). 

ft a bs,A i s the mass absorption coefficient, which we take 
from Draine (2003). The Draine (2003) model assumes 
that the dust consists of a mixture of carbonaceous grains 
and amorphous silicate grains, with dust grain size dis- 
tributions consistent with the wavelength-dependent ex- 
tinction in the local Milky Way, with By = 3.1 (Wein- 
gartner & Draine 2001). This model also assumes that 
most of the dust is heated by the diffuse radiation field, 
while the rest is heated by luminous stars with intense 
starlight. 

We compute the masses at A = 500/im, in order to 
minimize the dependence on the temperature, although 
the uncertainties of the flux densities are larger than 
at shorter wavelengths. At 500/im, the mass absorp- 
tion coefficient is K a b s = 0.95 cm 2 /g (which is « 20% 
lower than the value quoted in Li & Draine (2001)). Im- 
portantly, note that the absorption coefficients in these 
models have approximately k oc A -2 , while many of the 
galaxies have an emissivity closer to j3 » —1.5, which 
entails a wavelength dependence of the estimated dust 
masses. Masses estimated at A = 250/mi and 350/xm are 
lower by « 0.5 dex and 0.25 dex, respectively. 

For the galaxies lacking SPIRE fluxes, NGC 1404 and 
DDO 154, we estimate dust masses using the procedure 
outlined in Draine & Li (2007, Section 9.5), using the 
galaxies' fluxes at 8, 24, 70, and 160 /an. The Draine & 
Li (2007) model is an updated version of the one devel- 
oped by Li & Draine (2001). These dust masses are less 
reliable than the ones estimated from dust temperatures 
because the submm SEDs are less strongly constrained. 
The quoted errors of all of the dust masses consist only 
of the formal errors, and do not include systematic un- 
certainties. 

The dust masses of the KINGFISH galaxies fall be- 
tween 10 4 — 10 8 Mq. The dwarf and irregular galax- 
ies have the lowest masses, while early-type spirals such 
as NGC 1097 and NGC 7331 have the highest masses. 
Consistent with Masters et al. (2010a), luminous disk- 
dominated spirals such as NGC 4254, NGC 5457, and 
NGC 6946 have relatively large dust masses (Md us t ~ 
10 7 5 Mq). The dust masses and uncertainties are listed 
in Table [II 

Dust mass is strongly correlated with total infrared lu- 
minosity, which we analyzed in the previous section. The 
distribution of the dust /stellar flux ratio as a function of 
dust mass is qualitatively similar to Figure [31 but with 
larger uncertainties. 

Fifteen of our galaxies are included in the subsample of 
Draine et al. (2007) of galaxies that have SCUBA data. 
We have compared our dust masses to theirs, and our 
masses are systematically lower by ssO.2-0.4 dex. Some of 
their distances are shorter, and since distance appears in 



quadrature in Eqn. [3J accounting for this makes the dis- 
crepancy slightly larger. The discrepancy may be partly 
due to our single-temperature fit including the 70 fixa 
flux, which may have a contribution of small dust grains 
stochastically heated by starlight (Draine & Li 2007); 
galaxies with slightly larger temperature estimates would 
then have slightly smaller masses. The wavelength de- 
pendence of dust emissivity, mentioned above (and see 
Paradis et al. 2009), could also be an important factor. 
Two galaxies with large discrepancies are NGC 5713, 
for which Draine et al. (2007) obtained a mass nearly 
0.9 dex larger and for which a LMC dust model was fa- 
vored, and NGC 4631, for which they obtained a mass 
0.8 dex larger although the result depends on the allowed 
range of starlight intensity. Our dust mass for NGC 3077 
is slightly larger than that estimated by Walter et al. 
(2011), which we believe is mostly due to their smaller 
aperture size. A detailed analysis of the dust masses and 
other dust properties of KINGFISH galaxies will be pre- 
sented in Dale et al. (in prep.) and Aniano et al. (in 
prep.). 

4.2.2. Estimating M* 

We estimate the stellar masses from Zibetti et al. 
(2009), using optical and near-IR colors with if-band 
luminosity. In particular, they combine stellar popula- 
tion synthesis (SPS) models with simple prescriptions 
for dust attenuation. They use an updated version of the 
Bruzual & Chariot (2003) SPS models, which include re- 
vised prescriptions for the thermally pulsing asymptotic 
giant branch (TP-AGB) evolutionary phase, with a two 
component star formation history, consisting of a contin- 
uous, exponentially declining mode with random bursts 
superimposed. The Zibetti et al. (2009) model outputs 
stellar mass-to-light ratios logM*/ Lh(B — V, V — H). 
and for the LVL galaxies, which have ugriz-b&nd SDSS 
data, we use log Af* /£#((/ — i,i — H). In other words, 
the stellar masses are calculated with the following: 

logAf* = \ogL H + /(c pt,c optN i R ), (4) 

where c op t is the optical color (B — V or g — i) and c op tNm 
is the optical-NIR color (V — H or i — H). As a function 
of two colors, the model's mass-to-light ratios typically 
have 0.1-0.2 dex scatter. Before computing the masses, 
we have corrected the colors for foreground Galactic ex- 
tinction, using values taken from Schlegel et al. (1998). 
We assume a universal Kroupa (2001) initial mass func- 
tion (IMF). The assumed IMF affects the inferred masses 
and SFRs systematically (e.g., for a Salpeter IMF, they 
would be higher by a factor of 1.8 (Borch et al. 2006) and 
1.5 (Calzetti et al. 2010), respectively), while the relative 
distributions are robust. 

Most of the resulting stellar masses are in the range 
10 7 - 10 11 Af©. The SO galaxy NGC 1316, whose SED 
is shown in Figure [TJ has the largest stellar mass. Like 
the dust masses, the irregular galaxies have the lowest 
stellar masses; however, NGC 3077 is an exception: it 
is a relatively massive dwarf galaxy in the M81 group, 
whose starburst activity may have been triggered by in- 
teractions with its neighbors (Walter et al. 2002). Note 
that other models (e.g., Bell et al. 2003; Sani et al. 2011) 
yield relatively similar but not identical stellar masses for 
these galaxies; we refer the reader to Zibetti et al. (2009) 
for a comparison with and discussion of other models. 
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Fig. 4. — Dust mass (left) and stellar mass (right) versus dust/stellar flux ratio (Eqn.[TJ. Red circles, green squares, and blue triangles 
indicate galaxies with E and SO, Sa to Scd, and Sd and later-type morphologies, respectively. Open points indicate galaxies for which we 
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4.2.3. Results 

First, we show in Figure 0] the dust mass and stellar 
mass as a function of the dust/stellar flux ratio. The 
two figures are similar, because the galaxies in this sam- 
ple have such strongly correlated dust and stellar masses 
(unlike high-redshift submm galaxies, for example, which 
have a more scattered correlation; Santini et al. 2010). 

The galaxies are marked by their morphology classes 
in Figure H Four early-types (NGC 1404, NGC 584, 
NGC 1316, NGC 1291) and an earlier-type spiral 
(NGC 4594) are outliers, with low /dust//* and large dust 
and stellar masses. It is also interesting that /dust//* 
is correlated with the dust and stellar masses for the 
dwarf/irregular and late-type spiral galaxies, with Spear- 
man rank correlation coefficienttEH of r s — 0.52 and 0.51 
(which implies approximately 95% significance), respec- 
tively. Most dwarf galaxies lack substantial dust emis- 
sion (e.g., Walter et al. 2007), but the few of them with 
/dust//* > 0.3 (mostly Sd galaxies and NGC 3077) have 
relatively large dust and stellar masses. 

In both figures, but especially in the plot showing the 
stellar mass dependence, (Fig.|3b), there appears to be a 
transition between two populations of galaxies. In par- 
ticular, there appears to be a "transition" stellar mass at 
M* ~ 1O 1O M0, such that less massive galaxies follow a 
steep relation with dust / stellar flux, while more massive 
galaxies occur on a flat or slightly negative relation. We 
will later show that a similar transition appears to oc- 
cur between specific star formation rate and dust /stellar 
flux (see especially Fig. [TJj in Section T4.4p . This mass 
scale of 10 10 M Q is similar to the transition stellar mass 
determined by Kauffmann et al. (2003; see also Schimi- 

21 The Spearman rank correlation coefficient may have a value 
between —1 and 1. A positive (negative) value indicates an 
(anti)correlation, and a value of indicates no correlation. 



novich et al. 2007), above which galaxies have high stellar 
mass surface densities, high concentration indices typical 
of bulges, old stellar populations, and low SFRs and gas 
masses. Intriguingly, the maximum /dust//* occurs in 
the transition region, and perhaps NGC 1482 is a can- 
didate for a transition galaxy in this context; however, 
note that some other galaxies in this region are peculiar, 
and in some cases are interacting with a neighbor. In any 
case, it is interesting that the transition in /dust//* oc- 
curs at a similar stellar mass as the previously observed 
transition in D„(4000) (4000-Abreak strength, quantify- 
ing the star formation history), /i* (stellar mass surface 
density), and specific SFR. 

Next, we compare the ratio of dust /stellar mass to 
the dust /stellar flux ratios in Figure [5] Note that the 
Mdust/M* distribution is different than that of /dust//*: 
the /dust//* distribution is centrally peaked with small 
and large outliers (see Figure [2]), while the majority of 
the Mdust/M* ratios are between 10~ 4 and 10~ 2 , with 
the aforementioned early-types and NGC 4594 having 
lower values. NGC 584 in particular has a stellar mass 
that is more than five orders of magnitude larger than 
its dust mass. 

As stated in Table [21 the mean log dust/stellar mass 
ratio is — 2.95±0.09 for the sample, which is smaller than 
the low-redshift value given by Dunne et al. (2011), but 
the results appear to be consistent given the different 
sample selections. Dunne et al. find that the dust/stellar 
mass ratio depends on metal enrichment, which we also 
discuss below, and increases with increasing redshift (see 
also Bussmann et al. 2009). 

For the general galaxy population in KINGFISH, 
there is only a weak correlation between Mdust/-M* and 
/dust//* (Spearman rank r s — 0.30, less than 95% sig- 
nificance), with substantial scatter. Some of the scatter 
is likely due to systematic uncertainties in the dust and 
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stellar masses, as well as to the errors of the 500/im flux 
densities, which were used to estimate Mdust- Spiral and 
dwarf galaxies appear to populate the same locus on the 
plot. DDO 154 and DDO 165, both dwarf galaxies, are 
exceptions, with particularly small stellar masses. Three 
Sa galaxies — NGC 4594 in the lower left of the plot and 
NGC 1482 and NGC 2798 at the high /dust//* end— are 
also outside of the central locus of the plot, and appear 
to follow the trend of early-type galaxies. 

More clearly than the rest of the sample, the 
early-types exhibit a strong correlation between the 
dust/stellar ratios (Spearman rank r s = 0.56, 95% signif- 
icance). An interesting question to ask is why the early- 
types have such a large range of dust/stellar flux and 
mass ratios, spanning three orders of magnitude. One 
possibility is that these galaxies have a wide range of 
radiation field intensities and star formation rates (see 
Section 14.4)) , as well as widely varying gas / stellar ratios 
(Young et al. 2009). 

In particular, all but two of the early-types are lentic- 
ular galaxies, and many of them have similar masses and 
TIR luminosities. If we look at one pair of SOs, NGC 1316 
and NGC 1291 have extremely low dust/stellar ra- 
tios, and appear to be similar to the elliptical galaxies 
(NGC 1404 and NGC 584). Looking at another pair of 
SOs, NGC 1266 and NGC 1377 have particularly high 
dust/stellar ratios. For NGC 1266, we believe this could 
be due to an active nucleus (see Dale et al. 2007; Smith 
et al. 2007) which could be heating dust in the central 
region; however, this galaxy is likely an exception, as 
galaxies in the sample were selected not to have luminous 
AGN (Kennicutt et al. 2003). For NGC 1377, Roussel et 
al. (2006) argues that it is undergoing an opaque nascent 
burst of star formation, and the intense radiation field is 
significantly heating the dust. 

These two SO's and the Sa's mentioned previously, 
NGC 1482 and NGC 2798, are most offset from the trend 
of the rest of the sample. They are relatively bright in 
the far-IR, but this is not accompanied by a large dust 
mass. The offset may be explained by their relatively 
warm dust temperatures, which we discuss in Section r4.3l 
In principle, large stellar mass-to-light ratios could also 
contribute to offsets in this direction, but the M*/L of 
these galaxies are not particularly large. 

We emphasize that the outliers and substantial scat- 
ter evident in Figure [5| highlight the danger of using 
/dust//* as a proxy for Mdust/M*. The two quantities 
are certainly related, but they probe different physical 
processes, with different dependencies on a galaxy's star 
formation history and history of dust production and de- 
struction. Considering the selection criteria of the KING- 
FISH survey, only one fourth of which is composed of 
faint dwarf galaxies, the "true" scatter between these 
dust/stellar ratios is probably even larger. For the rest 
of this paper, we continue to focus on relations between 
/dust//* and other galaxy properties, in order to keep 
the analysis as empirical as possible; subsequent KING- 
FISH papers will provide more detailed analyses of the 
galaxies' dust and stellar masses. 

We also briefly note that the dust/stellar mass ra- 
tios vary slightly with bar strength: the mean ratio for 
strongly barred galaxies (SAB or SB) is log Md us t /M* = 
—2.84 ± 0.58 (rms), while galaxies with weak or no bars 



have a slightly lower mean ratio, —3.06 ± 0.74 (rms). 
This dependence may be related to the fact that strongly 
barred galaxies have warmer dust temperatures in their 
inner regions, likely due to bar-induced star formation 
in these regions (Engelbracht et al. 2010). Nonethe- 
less, the variation of Md ust /M» among barred and non- 
barred galaxies is much greater than the difference be- 
tween them. 

The metallicity dependence of the dust/stellar corre- 
lation (Figure \5jp) is similar to the morphology depen- 
dence, because of the relation between morphology and 
metallicity (e.g., Moustakas et al. 2010). Nonetheless, 
there are a few differences. The correlation is strongest 
for metal-rich galaxies (r s — 0.80, 99.9% significance), 
while some galaxies with intermediate metallicities have 
lower dust/stellar mass ratios than one might expect, 
given their flux ratios. These are SO and Sa galaxies 
with larger than average star formation rates (see Sec- 
tion I4.4[) ; their new stars could be significantly heating 
the dust and increasing the FIR emission, without sig- 
nificantly increasing the dust mass. 

In addition, the fact that the dust/stellar flux ra- 
tio is more metallicity-dependent than the dust/stellar 
mass ratio (comparing the metal-poor and intermediate- 
metallicity galaxies) suggests that the scatter between 
them is partly due to metallicity. This appears to be 
consistent with Groves et al. (2008), who found in their 
SED models that the IR peak shifts to longer wavelengths 
and becomes broader with increasing metallicity, due to 
the increasing dust column and the increasing mechani- 
cal luminosity of starburst regions. The residuals are also 
significantly correlated with dust temperature, suggest- 
ing that the scatter is also partly due to temperature, and 
the geometrical distributions of dust and stars within the 
galaxies. This is not surprising: dust mass is inversely 
related to dust temperature (Eqn. [3]) , and as we show in 
the next section, the dust/stellar flux ratio is positively 
correlated with temperature, which implies that galaxies 
with warmer dust tend to have more dust emission, and 
tend to be located towards the lower right of Figure [5] 

Lastly, we note that the metallicity dependence of the 
dust/stellar mass ratio here is qualitatively similar to 
that of the dust/gas mass ratio in Draine et al. (2007), 
who studied a similar sample of galaxies. For the galax- 
ies with submm data, Draine et al. (2007) find that the 
dust /gas ratio gradually increases by up to 0.3 dex over 
the range of gas-phase metallicity. We also find that 
our metal-poor galaxies tend to have lower dust/stellar 
flux ratios than intermediate-metallicity galaxies. Nev- 
ertheless, there is no significant metallicity dependence 
of the dust/stellar mass ratios. Therefore, the scatter 
in Figure [5] may be partly explained by the differential 
dependence on metallicity. 

Furthermore, since lower-metallicity galaxies have 
lower dust/gas mass ratios but not necessarily lower 
dust/stellar mass ratios, it is possible that they have 
lower stellar mass for a given amount of gas — that is, 
that their star formation is less efficient than more metal- 
rich galaxies. This is consistent with Lee et al. (2006), 
who found that metallicity is inversely related to the 
gas/stellar mass ratio. Note, however, that specifically 
for massive galaxies, Schiminovich et al. (2010) argue 
that star formation efficiency is independent of stellar 
mass and stellar mass surface density. 
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Fig. 5. — Dust/stellar mass ratio versus dust/stellar flux ratio. Left plot: red circles, green squares, and blue triangles indicate galaxies 
with E and SO, Sa to Scd, and Sd and later-type morphologies, respectively. Open points indicate galaxies for which we do not have 
detected SPIRE fluxes and dust temperatures; these dust masses may have additional systematic uncertainties. Right plot: red circles, 
green squares, and blue triangles have the highest, intermediate, and lowest oxygen abundances, from Moustakas et al. (2010); open points 
indicate abundances estimated from the luminosity-metallicity relation. 



4.3. Far-Infrared Colors 

We now examine the far-infrared colors of the galaxies 
in our sample, using MIPS and SPIRE bands. 'Warmer' 
far-IR colors are often associated with higher tempera- 
tures of small or large dust grain components, depend- 
ing on the wavelengths (e.g., Li & Draine 2001; Boselli 
et al. 2010b). Low / 1/ (70/im)// ! ,(160/im) color, for ex- 
ample, indicates that the far-IR SED peaks at a long 
wavelength, likely due to a large fraction of cold dust. 

In Figure [SJ we show the correlation between the 
dust/stellar flux ratio and /„(70/im)// 1/ (160/im). We 
also show the cor relat ion with dust temperature, de- 
scribed in Section 14.21 which uses these bands as well 
as the SPIRE flux densities. The correlations between 
the dust /stellar flux ratio and SPIRE colors (not shown) 
have similar trends, but with more noise. As stated in 
Section |4. 2. 1[ the temperatures depend on the assumed 
emissivity (e.g., Bendo et al. 2003), though the rela- 
tive distribution of temperatures is robust. NGC 1377 
is the galaxy with the warmest dust temperature, as 
mentioned in the previous section, as well as the largest 
/„(70/im)// I ,(160/im) color. The trends in both figures 
are similar because an increased /„(70yLtm)//„(160/im) 
color generally corresponds to a warmer dust tempera- 
ture, even though we are simply fitting a single temper- 
ature modified blackbody. 

In both figures, there is a correlation, but with consid- 
erable scatter (Spearman rank r s = 0.33 for the color plot 
and r s = 0.36 for the Td us t plot, with 95% significance). 
In particular, only early-type galaxies (r s = 0.68 and 
r s = 0.78 for the two plots) and spirals (r s = 0.66 and 
r s = 0.54) exhibit a significant (3-<r) correlation: galax- 
ies with relatively high dust emission also have relatively 
high far-IR colors and warm temperatures. 



In addition, at a given /dust//* ratio, early- types have 
warmer colors and dust temperatures (by up to 5 K) than 
spirals. (The single exceptional early-type is NGC 1291, 
which Hinz et al. (in prep.) show has an outer ring that 
is dominated by much cooler dust than the rest of the 
galaxy.) Perhaps this occurs because the early- types in 
KINGFISH tend to have slightly smaller dust masses, 
which require less of a heating source (e.g., stars) to 
significantly raise their temperatures. It has long been 
known that dust grains may be heated by old stars (see, 
e.g., Helou 1986; Sauvage & Thuan 1994; Calzetti et al. 
1995; Kennicutt 1998; Draine & Li 2001); however, the 
clearly different Td us t — /dust//* trend for early- types 
and spirals is a new result. The SO's with surprisingly 
large /dust//* also have relatively warm dust tempera- 
tures, as well as slightly larger specific star formation 
rates, which we show in Section 14.41 It appears that 
some SO's in the sample (as well as a couple Sa's) are un- 
dergoing a period of obscured star formation, and may 
be similar to the early-types examined by Shapiro et al. 
(2010) and Wei et al. (2010), which have molecular gas 
possibly accumulated through minor mergers, as well as 
slightly elevated star formation efficiencies. Nevertheless, 
most of the early- types with /dust//* < 1 have very lit- 
tle ongoing star formation and lower stellar fluxes than 
their spiral counterparts, and yet they still have slightly 
warmer temperatures. Hence, these early-types may sim- 
ply have more intense radiation fields. Early- types have 
many stars in their bulge components (as indicated, for 
example, by their 3.6 /im surface brightnesses), and these 
regions tend to have warmer dust temperatures (Engel- 
bracht et al. 2010) in spite of the lower star formation 
rates, so we argue that a more intense interstellar ra- 
diation field is the most likely explanation of the differ- 
ent trends for late- and early-types. This issue is cur- 
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Fig. 6. — Left plot shows MIPS /i/(70/Ltm)//i/(160/im), a proxy for dust temperature, versus dust/stellar flux ratio. Right plot shows the 
dust temperature itself, estimated from a fit to the far-IR SED with a single-temperature modified blackbody, vs. dust/stellar flux ratio. 
Galaxies with E and SO, Sa to Scd, and Sd and later-type morphologies are indicated by red circles, green squares, and blue triangles, 
respectively. 



rently being investigated further, using models of the 
dust masses and starlight intensities (Aniano et al., in 
prep.). 

Finally, the dwarf and irregular galaxies are scat- 
tered in the figures, although as noted by Walter et 
al. (2007) and Dale et al. (2007), some dwarfs have 
high / J/ (70/xm)// J ,(160/xm) ratios, indicating strong over- 
all heating of the dust grain population. 

4.4. Star Formation Rate 

4.4.1. Estimating SFRHa+iA^m and SFRf\jv+tir 

Finally, we turn to the star formation rates (SFRs) of 
the galaxies in our sample. One way to estimate SFRs 
is to combine Ha and mid-IR (specifically, 24 fjxa) lumi- 
nosities (Calzetti et al. 2007; Zhu et al. 2008; Kennicutt 
et al. 2009). We use the calibration proposed by Calzetti 
et al. (2010), with scatter < 0.2 dex, for galaxies with a 
wide range of metallicities: 

SFR (M /yr) = C Hq [L(Ra) ohs + aiL(24)] (5) 
if £(24) < 4 x 10 42 erg/s, 
= C Hct [£(Ha) obs + a 2 £(24)] (6) 
if 4 x 10 42 < £(24) < 5 x 10 43 erg/s 

where the luminosities are in units of ergs -1 , Ch q = 
5.45 x lO -42 (M yr -1 )/(ergs -1 ), ai = 0.020, and a 2 = 
0.031. The first of these equations (J3J) is calibrated for 
normal galaxies (Kennicutt et al. 2009), while the latter 
([6]) is calibrated for HII regions and starbursts (Calzetti 
et al. 2007). Based on the 24 /im luminosity, most of the 
KINGFISH galaxies (42/62, and 41/54 with Ha fluxes) 
qualify as 'normal' galaxies. 

For comparison, we also estimate SFRs using IR and 
UV luminosities, which is a complementary way to ac- 
count for both obscured and unobscured star formation 



(e.g., Zheng et al. 2007; Buat et al. 2007). We use the 
SFR parameterization of Hao et al. (2011), which com- 
bines far-UV (1500 A) and TIR (Eqn. [2} luminosities, 
and which updates the calibration in Dale et al. (2007): 

SFR{Mq/ Y t) = C FlJV L fuv (1 + a L Tm /L FlJV ), (7) 

where the luminosities are in units of ergs -1 , Cfuv = 
4.30 x lO -44 (M yr -1 )/(ergs -1 ), and a = 0.475. Like 
the stellar masses, we have assumed a Kroupa (2001) 
IMF for the SFRs. 

There is some dependence on the assumed star for- 
mation history (SFH), such as a SFH with a short re- 
cent starburst or with constant or declining star forma- 
tion over a long time-scale, but the largest systematic 
uncertainty is the assumed (universal) IMF. There also 
are other systematic uncertainties due to metallicity and 
AGN activity and to assumptions about the attenua- 
tion correction (see Kennicutt et al. 2009; Calzetti et 
al. 2010). 

Both star formation rates for the KINGFISH sample 
are listed in Table [TJ except for the galaxies lacking Ha 
or far-UV fluxes. The errors listed in the table are only 
the formal errors due to the flux densities; systematic 
uncertainties contribute at least 0.2 dex of additional un- 
certainty. The SFRs of the galaxies in the sample range 
from 10 -3 to 10 M yr -1 . 

The two SFR estimates are generally consistent, within 
0.3 dex of each other. There are a few galaxies with larger 
discrepancies, such as NGC 1377 (although we had to 
perform an uncertain extrapolation for the far-UV flux) , 
NGC 3351, NGC 4736, and NGC 6946 (although it has 
substantial diffuse background emission). IC 342, which 
is close to the Galactic plane (I = 138°, b = 10°), is dom- 
inated by obscured star formation, and the discrepancy 
between its SFRs appears to be due to a relatively low 
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24 /j,m flux. Consequently, its FUV+TIR SFR is likely 
more accurate. There are also discrepancies for low-mass 
dwarf irregular galaxies with very low SFRs (IC 2574, 
DDO 154, and DDO 165; see Walter et al. 2007), which 
are difficult to determine accurately. 

Our SFR and stellar mass estimates (discussed in Sec- 
tion I4.2.2[) can be compared to those obtained from 
model fits to the SEDs. For example, Johnson et al. 
(2007) used the stellar population models of Bruzual & 
Chariot (2003) and the dust models of Witt & Gordon 
(2000). Using similar models, we generally find excellent 
agreement for the stellar masses, and good agreement 
for the SFRs, which are consistent within 0.3 dex. Two 
exceptions are NGC 3351, for which the SED fit yields 
a SFR that is 0.6 dex larger than our estimate using 
Ha + 24/im luminosities, and NGC 3521, for which the 
fit yields a SFR that is 0.6 dex smaller. 

We also compare to da Cunha et al. (2008), who use 
the same stellar population models with exponentially 
declining star formation histories with bursts, and with 
the dust emission modeled as the sum of modified black- 
bodies at different temperatures. We also find excel- 
lent agreement for the stellar masses, but the SFRs of 
nine galaxies (out of 57) are statistically inconsistent 
and deviate by more than 0.4 dex. The most extreme 
case is NGC 3190, for which da Cunha et al. obtain 
log SFR = —1.23, more than 0.8 dex lower than our 
estimates. Finally, Noll et al. (2009) also estimated stel- 
lar masses and SFRs of a sample of galaxies, 32 of which 
are in our KINGFISH sample. They used the Maraston 
(2005) stellar population models, assuming either con- 
stant star formation or an exponential decreasing SFR. 
Their SFRs are mostly consistent with ours, except for 
NGC 4536 and NGC 4736, for which they obtain values 
much larger than our SFR(ffa + 24/zm) but similar to 
our SFR(TIR+FUV). We conclude that our SFRs are 
generally reliable, although for some galaxies they are 
difficult to determine accurately within a factor of two. 

4.4.2. Results 

We now show the correlation between specific SFR, or 
SFR per unit stellar mass, and dust/stellar flux ratio in 
Figure [3 using both SFR estimators. The advantage of 
using SFR/M* over SFR is that it allows us to fairly 
compare galaxies with a wide range of stellar mass. 

In general, there is a weak trend such that specific SFR 
increases with the dust/stellar flux ratio. In other words, 
galaxies with relatively bright dust emission also tend to 
have more star formation. The scatter appears to be 
larger using SFR// Q+ 24 Mm , although a statistical anal- 
ysis indicates that its correlation is of similar strength 
(Spearman rank r a = 0.43, versus 0.30 for SFRfuv+tir)- 
Early-type and spiral galaxies exhibit particularly tight 
correlations (r s w 0.7 and r s w 0.65, respectively, with 
99% significance). The scatter in the figure may be due 
to a number of factors, such as different contributions 
from obscured and unobscured tracers of star formation. 

The correlation between the dust/stellar flux ratio and 
specific SFR is interesting in the context of its corre- 
lation with stellar mass and metallicity (see Figs. 0Jd 
and [SJd), considering that some authors have recently 
argued for a 'fundamental plane' relating stellar mass, 
SFR, and metallicity of emission-line galaxies (Mannucci 
et al. 2010; Lara-Lopez et al. 2010). Many of the KING- 



FISH galaxies would lie on this plane, but with substan- 
tial scatter and with outliers, among which would be 
the early- type galaxies lacking prominent emission lines. 
These three galaxy properties — stellar mass, SFR, and 
metallicity — are certainly related to the dust/stellar flux 
ratio, but the inter-relations between star formation, stel- 
lar mass growth, dust production and metal enrichment 
are complex. 

We can compare our results in Figure [7] to da Cunha et 
al. (2010), who studied similar properties of SDSS galax- 
ies by modeling their SEDs. They similarly obtain a cor- 
relation between Md UB t/M t and SFR/Af,, although their 
correlation appears to be slightly stronger. They argue 
that stellar mass is not the main driver of this correlation. 
Nonetheless, we find that plotting the dust /stellar flux 
ratio against SFR (without normalizing by M*) yields 
a shallower correlation with more scatter, da Cunha et 
al. (2010) also find that SFR/M* is strongly correlated 
with the dust-to-gas ratio and the fraction of -Ltir con- 
tributed by dust in the ambient ISM. In any case, our 
results appear to be consistent with theirs, with small 
differences likely due to sample selection (their sample is 
dominated by star-forming galaxies) and to the fact that 
their SEDs do not probe wavelengths longer than 100/im, 
which could yield underestimates of the dust fluxes and 
masses for some galaxies. 

The dwarf and irregular galaxies in Figure [7] de- 
part from the positive correlation between specific SFR 
and dust/stellar flux ratio that we see for earlier- type 
galaxies (r s = —0.33 for SFR/f Q+ 24 Aim and —0.70 for 
SFRfuv+tir; the latter anti-correlation is statistically 
significant). Somewhat similarly, Dale et al. (2007) find 
that dwarf/irregular galaxies have specific SFRs that de- 
crease with the infrared-to-ultraviolet ratio (their Fig. 
10). The dwarf galaxies' location in Figure [7] is most 
likely due to their small stellar masses, and their limited 
dust content (Walter et al. 2007); only a small fraction 
of dwarfs have large SFRs (Lee et al. 2009) . 

Late-type and early-type galaxies have a similar rela- 
tion between SFR/M* and /dust//*- Nonetheless, it is 
interesting that, according to Figure [H late- types have 
cooler dust temperatures (by up to 5 K) at a given 
dust/stellar flux ratio. For example, some SO and Sa 
galaxies have similar specific SFRs and /dust//*, but the 
SOs have slightly warmer dust and lower dust masses. 
Perhaps in some of these galaxies, a small amount of star 
formation can more efficiently heat the dust; it is more 
likely, however, that the heating by the general ISRF, 
as well as by massive stars, is sufficient to heat the dust 
grains slightly more (e.g., see Draine & Li 2007). 

It is also interesting that approximately half of the 
early-type galaxies in the KINGFISH sample still have 
on-going star formation and dust heating, and are more 
similar to some of their late-type counterparts than to 
the passive early-types. This is contrary to the view that 
galaxies with large bulge components quickly 'quench' 
their star formation. Fabello et al. (2011) argue against 
this view, showing that early-types with large bulge-to- 
disk ratios do not have lower atomic gas contents or 
star formation efficiencies than similar later type galax- 
ies. Crocker et al. (2011) account for both atomic and 
molecular gas and obtain a similar result, such that 
star-forming early-type galaxies lie in a similar range of 
the Schmidt-Kennicutt relation as normal star-forming 
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Fig. 7. — Correlation between specific star formation rate and dust/stellar flux ratio. The red circles, green squares, and blue triangles 
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Section [3}. 



galaxies. Bulge components may be a necessary con- 
dition for suppressing star formation, but they are not 
sufficient (Bell 2008). 

A number of KINGFISH galaxies appear to be similar 
to 'passive disks' and 'red spirals^! recently studied by 
other authors (e.g., Wolf et al. 2009; Skibba et al. 2009; 
van der Wei et al. 2009; Masters et al. 2010b; Bundy et 
al. 2010). Ten galaxies in the sample have dust masses 
greater than 10 8 Af Q , and all of them are spirals. The 
majority (7/10) of the galaxies with /dust//* > 0.75 are 
spirals as well. Of the spiral galaxies with large dust 
masses or large dust/stellar ratios, some have low spe- 
cific SFRs, such as the earlier-type spirals NGC 4725 
and NGC 1512. These may be examples of 'transition' 
galaxies discussed by Masters et al. (2010b) and Bundy et 
al. (2010), galaxies between star-forming disk-dominated 
galaxies and passive early-types. 

Nonetheless, large dust mass is not a sufficient con- 
dition for declining SFRs, because there are dusty star- 
forming galaxies such as the late-types NGC 4254 and 
NGC 4631. Earlier- type spirals with obscured star for- 
mation, such as NGC 1097, NGC 2146, NGC 2798, and 
NGC 7331, have also been described as transition galax- 
ies (Zhu et al. 2011). In addition, a galaxy need not 
have a dominant bulge to have passive star formation: 
NGC 5055 and NGC 4321 are dusty galaxies with very 
low specific SFRs, but with Sbc morphologies. Finally, 
the KINGFISH sample also contains massive galaxies 
with 'pseudobulges' (i.e., lacking classical bulges), like 
NGC 5457 and NGC 6946, which are difficult to explain 
with current models of galaxy formation, in which most 

22 'Red spirals' are spiral galaxies on the red sequence, as defined 
by optical colors. They include both passive spirals as well as 
spirals with obscured star formation (Masters et al. 2010b). 



massive galaxies experience major mergers (Kormendy et 
al. 2010); however, they both have a very large gas sup- 
ply (Walter et al. 2008) with which they could potentially 
form stars, and NGC 6946 may also have a particularly 
high star formation efficiency in its spiral arms (Foyle 
et al. 2010). It appears that, while some galaxies may 
be undergoing a transition to passive star formation and 
earlier- type morphologies, the transition is blurry and 
some galaxies evolve differently than others. 

5. EXTRAGALACTIC BACKGROUND LIGHT 

Our dust/stellar flux ratios quantify the relative con- 
tributions of emission from dust and stars of individual 
galaxies. The extragalactic background light (EBL) is 
a related concept, because it can be used to quantify 
the contributions of stellar- and dust-dominated galax- 
ies to the optical and infrared background. Rather than 
examining SEDs within galaxies, analyses of the EBL in- 
volve integrating the light from many background galax- 
ies. Studies of the EBL distinguish between optical light 
from stars and IR light from dust, so it is useful to com- 
pare these studies to our analysis of the stellar and dust 
emission of individual galaxies, although because of the 
small size and incompleteness of the KINGFISH sam- 
ple, we cannot make strong conclusions based on these 
comparisons. 

A number of authors have investigated the EBL and 
have quantified it at different wavelengths (e.g., Hauser 
et al. 1998; Hauser & Dwek 2001; Dole et al. 2006; 
Bethermin et al. 2010; Kneiske & Dole 2010). A re- 
cent study (Dommguez et al. 2011) has attempted to 
distinguish between the contributions of different types 
of galaxies to the EBL. There has also been some theo- 
retical work on explaining the intensity of the EBL, in 
the context of the evolving stellar mass density and SFR 
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density (e.g., Chary & Elbaz 2001; Primack et al. 2005; 
Nagamine et al. 2006; Fardal et al. 2007). For a recent 
comparison of some models and constraints on the IR 
EBL, see Orr et al. (2011). The EBL is the dominant 
radiant energy in the universe aside from the cosmic mi- 
crowave background, and it is thought to be mostly due 
to stars, dust, and AGN in galaxies (e.g., Hauser & Dwek 
2001), while faint diffuse emissions can represent only a 
small fraction of the integrated energy (Dole et al. 2006). 

Some have attempted to compare the cosmic infrared 
background (CIB) and the cosmic optical background 
(COB). Hauser & Dwek (2001) found that « 52% of the 
EBL is contributed by direct starlight (and the rest ab- 
sorbed and re-emitted by dust at A > 3.5/im), but with 
large uncertainties, while Dole et al. (2006) found that 
the ratio COB/CIB ranges from 0.7 to 1.5 (making the 
demarcation at A = 8/um). Perhaps not coincidentally, 
some of the galaxies in the KINGFISH sample have sim- 
ilar ratios of /dust//*, which is analogous to CIB/COB. 

The galaxies that dominate the COB and CIB may 
constitute different populations; for example, Dole et al. 
(2006) argue that the CIB is mainly composed of lu- 
minous infrared galaxies (LIRGs) at z ~ 1. Galaxies in 
the KINGFISH sample with bright TIR luminosity, large 
masses and specific SFRs, such as NGC 1482, NGC 2146, 
NGC 2798, and NGC 7331, could be considered to be 
examples of low-redshift counterparts of these galaxies, 
although they are not as luminous as LIRGs. KINGFISH 
also includes galaxies which may be counterparts of ob- 
jects dominating the optical EBL, such as NGC 3521 and 
NGC 5055, which are massive galaxies with significant 
star formation. 

Lastly, note that the mean stellar/dust flux ratio of the 

galaxies in our sample is /dust//* = 0.55, and split by 
morphological type, the means are 0.28, 0.69, and 0.54 
for dwarfs, spirals, and early-types, respectively. The 
summed /dust//* ratio can be more directly compared 
to the CIB/COB estimates of Hauser & Dwek (2001) 
and Dole et al. (2006). This quantity, which is domi- 
nated by the more luminous galaxies, is £(/d U st)/£(/*) = 
0.44 ± 0.01 for the whole sample. This can also be ex- 
pressed as £(/dust)/£(./* +/dust) = 0.30 ±0.01. Split by 
type, this dust fraction is 0.34, 0.32, and 0.07 for dwarfs, 
spirals, and early-types, respectively. KINGFISH does 
not contain a statistically complete sample; in general, 
dwarf galaxies dominate by number, although they are 
typically faint (e.g., Lee et al. 2009). Many early-types 
are luminous and stellar-dominated, but they are rela- 
tively few in number (e.g., Dale et al. 2009). Our re- 
sults appear to indicate that, at least for nearby galax- 
ies, spirals may constitute the most important contri- 
bution, for which their stars contribute ~ 68% of the 
light and their dust contribute ~ 32%. This is approxi- 
mately consistent with Soifer & Neugebaucr (1991), who 
analyzed a flux density-limited sample of local galax- 
ies with IRAS and obtained a total dust/stellar frac- 
tion of « 23%, which is smaller than our result, prob- 
ably due to their limited coverage in the far-IR. In any 
case, the local fraction is smaller than the high-redshift 
infrared EBL fraction estimated by Dole et al. (2006; 
CIB/(COB+CIB)« 40 - 60%), which indicates that the 
IR output of galaxies evolves with time and was larger 
in the past, consistent with studies of the evolution of IR 



luminosity functions (e.g., Le Floc'h et al. 2005; Chary & 
Pope 2010; Murphy et al. 2011a; cf., models in Fontanot 
& Somerville 2010). Nevertheless, a larger and more 
complete low-redshift catalog would be required to in- 
vestigate this issue further. 

6. SUMMARY 
We now summarize our main conclusions: 

• The dust/stellar flux ratios estimated from global 
galaxy SEDs are correlated with total-infrared lu- 
minosity, in a morphology and metallicity depen- 
dent way. Metal-poor dwarf galaxies tend to have 
faint IR luminosities, while spirals tend to have 
lower mctallicities and higher dust /stellar flux ra- 
tios than early-types. 

• Dust/stellar flux and dust/stellar mass ratios are 
correlated, especially for early-types. Late-types 
and dwarf galaxies show considerable scatter, 
partly due to the effect of metallicity on the flux 
ratio. Some of the scatter is also due to the depen- 
dence on dust temperature. 

• Most galaxies exhibit a trend such that, those with 
large dust /stellar flux ratios have warmer FIR col- 
ors and dust temperatures. Late-types tend to have 
slightly cooler temperatures (by up to 5 K) than 
early- types at a given dust/stellar flux ratio, while 
dwarf and irregular galaxies have more scattered 
temperatures. 

• We find that late- and early-type galaxies have spe- 
cific SFRs that are correlated with dust /stellar flux 
ratios: galaxies with more dust emission also tend 
to have more star formation. Combined with our 
previous result, we interpret this as evidence that 
ongoing star formation is sufficient to heat some 
of the dust in these galaxies, while other galaxies 
have more intense radiation fields, where the older 
stellar population likely contributes significantly to 
the dust heating. 

• The KINGFISH sample contains a number of dusty 
star-forming spiral galaxies as well as some passive 
spirals, whose limited star formation resembles that 
of some early- types. 

Finally, we note that our results could contribute con- 
straints for galaxy formation models, such as on the 
amount of dust production, metal enrichment, and star 
formation in different types of low-redshift galaxies. In 
addition, our results could be useful as a local benchmark 
for comparisons with high-redshift studies, such as stud- 
ies of submillimeter galaxies (e.g., Santini et al. 2010) 
and of the evolution of the dust and stellar content of 
galaxies (e.g., Dunne et al. 2011). 
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